Search for E 2ff phonon modes in MgB 2 single crystals by point-contact spectroscopy 
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The electron-phonon interaction in magnesium diboride MgB2 single crystals is investigated by 
point-contact (PC) spectroscopy. For the first time the electron coupling with E2 9 phonon modes is 
resolved in the PC spectra. The correlation between intensity of the extremely broad E2 3 modes in 
the PC spectra and value of the superconducting gap is established. Our observations favor current 
theoretical models for electron-phonon mediated superconductivity in MgB2 and they better match 
the harmonic phonons model. 

PACS numbers: 74.25.Fy, 74.80.Fp, 73.40.Jn 



Introduction. During the two year rush in study of su- 
perconductivity in MgB2 the electron-phonon coupling 
has been debated as a main source of the Cooper pair- 
ing. Many of theoretical calculations report similar find- 
ings: the Ei2g phonon modes, corresponding to the in- 
plane distortions of the boron hexagons, have an enor- 
mous coupling to electrons and are responsible for high 
T c |, | | |, |, UJ. Correspondingly, the main fea- 
ture of the calculated electron-phonon interaction (EPI) 
function a 2 F(uj) is a dominant maximum at the energy 
of stretching E 2ff modes within 60-70 mcV §], |, §, @, ff . 
However, in spite of the great progress in computation of 
a 2 F(uj), which is calculated even for various parts of the 
Fermi surface [Q, |), direct experimental verification for 
the a 2 F(uj) behavior is still lacking. 

The importance of the Fi2 g modes is evidenced by the 
Raman spectroscopy studies, which demonstrate strongly 
damped and broad optical modes with the E2 ff symmetry 
H ^, [HJ O. Additionally, an anomalous large broad- 
ening of the Fi2 g phonon modes along the T-A direction 
is recently observed in single crystals by inelastic X-ray 
scattering It is shown that the dominant contribu- 
tion to the linewidth is due to a strong electron-phonon 
interaction (EPI). But both the foregoing experiments 
deal with phonons and not directly with the a 2 F(uj) func- 
tion. 

The two well known fundamental experimental meth- 
ods for determining a 2 F(u) are point-contact (PC) and 
tunneling spectroscopy. Both of them were successfully 
utilized for observation and study of superconducting 
gap(s) (see, e. g., Refs. i@§@§i§). Nev- 
ertheless, only a few papers report the search for the 
phonon structure and EPI in PC [EOl El| or in tunneling 



spectra |22 . 

Bobrov et al. p0[ | communicated reproducible max- 
ima, including a maximum around 60 mV, in the PC 
spectra for different contacts prepared on the c-axis ori- 
ented MgB2 films. In the forthcoming paper Yanson et 



al. |2l| stressed that in the superconducting state the 
phonon structure observed in the PC spectra was due 
to the energy dependence of the superconducting order 
parameter [^3], For strong coupled superconductors 
this, so called, elastic term |24| exceeds the inelastic EPI 
contribution to the PC spectra determined directly by 



Consequently, to separate or to distinguish 



the elastic and inelastic effects measurements of PC spec- 
tra of MgB2 in the normal state are required. However, 
the resolution in the point-contact spectroscopy (PCS) is 
limited by temperature and at T > T c ~ 40 K it is above 
20 mV. On the other hand, the MgB2 thin films have a 
large critical field above 25 T [^6| to drive them in the 
normal state at helium temperature. 

This paper reports the investigation of MgB2 single 
crystals by PCS. The advantages of the single crystal as 
compared to the films in the context of the foregoing dis- 
cussion are a few times lower critical field, permitting su- 
perconductivity to be almost suppress by moderate field 
along the c-direction , the best sample quality (lower 
residual resistance or larger mean free path) and the pos- 
sibility of anisotropy study. 

The goal of the investigations is to recover the EPI 
function for MgB2 by means of PCS in order to elucidate 
the mentioned above issue about role of E2 9 modes in 
superconductivity of this compound. 

Experimental details. The single crystals of MgB2 were 
grown in a quasi-ternary Mg-MgB2-BN system at a pres- 
sure of 4 - 6 GPa and temperature 1400-1700 °C for 5 to 
60 min. In optimal conditions shiny yellow-colored single 
crystals with the large side of 0.2 - 0.7 mm were yielded. 
The details of the sample preparation are given elsewhere 
p7j . The estimated resistivity at 40 K is about 1 fiCt cm 
and the residual resistivity ratio RRR=5± 0.1. A sharp 
superconducting transition in resistivity is around 38.1- 
38.3 K with <5T c (10-90)% = 0.2 - 0.3 K. 

The PC's were established in situ at low temperatures 
by touching of a Cu or Au electrode by tiny piece (scrap) 
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of the Mgl$2 single crystal. That is because of the small 
size and the irregular shape the MgB2 crystal was used 
as a "needle". We always try to align the c axis of the 
MgB 2 single crystal with an applied magnetic field and 
the contact axis was along the largest side, presumably, 
perpendicular to the c axis in most cases. Due to lack of 
special micro-mechanics the accuracy of orientation was 
held by eye. No more control of the sample orientation 
was possible during the measurements. Therefore, infor- 
mation about the contact axis orientation with respect 
to the main crystallographic direction was gained from 
the analysis of the gap features. 

Results and discussion. Study of nonlinear conductiv- 
ity of metallic PC's allows the EPI function a 2 F(ui) to 
be recovered directly by measuring the second derivative 
of the I — V characteristic (see, e. g., Ref. EE 
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AR 1/2 a 2 PC (e)F(e)U =cV , 



(1) 



where A is the constant, R — dV/dl, Rq is the contact 
resistance at zero bias, ctp C F(uj) is the PC EPI func- 
tion, which differs from the Eliashberg EPI function by 
the presence of factor K = 1/2(1 — #/tan#), where 6 is 
the angle between initial and final momenta of scattered 
electrons [for the transport and Eliashberg EPI functions 
the corresponding factors are: K = (1 — cos#) and K=l, 
respectively]. From Eq.(Q) ap C (e) F(e) can be expressed 
via the measured rms signal of the first V\ ocdV/dl and 
the second V2 ocd 2 V/dI 2 harmonics of a small alternating 
voltage superimposed on the ramped dc voltage V: 
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It should be emphasized that to discard the influence of 
the superconducting peculiarities we have analyzed the 
PC spectra measured at a maximal magnetic field (as 
high as 9T) which is close to the upper critical field in 
the MgB 2 single crystal |§7|) and/or above T c <40K. 
The common features for all the spectra measured at low 
temperatures are the residual superconducting structure 
(the sharp maximum near zero bias) and the continu- 
ously raising structureless background above 25-30 mV. 
For a number of curves a smoothed maximum (rather a 
hump) between 40-80 mV was resolved. The examples of 
the measured dependences with the maximum are shown 
in Figs. 1-3. First of all we should stress that the struc- 
ture in the spectra shown above 30 mV is robust with 
varying the magnetic field in the range 6-9 T, while the 
superconducting structure below 30 mV depresses dras- 
tically by a magnetic field (see Fig. 1). Secondly, above 
30 mV the d 2 V/dI 2 (V) curves are reproducible for both 
polarities of the applied voltage. This is in support of the 
phonon nature of the maxima(um) in d 2 V / dl 2 (V) . From 
the minima position in the dV /dl curves in the insets of 
Figs. 1 and 2 measured at zero field and low temperature 
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FIG. 1: PC spectra for the MgB 2 -Cu contact (R = 7.2 fi) 
at T = 4.2K in a magnetic field 8T (dotted curve) and 9T 
(solid curve). The dashed curve is for the opposite polarity at 
9 T. Above 30 mV the curves practically coincide. The dash- 
dotted curve is recorded at zero field at 40 K. The modulation 
signal Vi(0) is 3mV. Inset: the dV / dl (V) curves for the same 
contact. The bottom two curves are for zero and 1.2 T field 
and 4.2K. 

we can estimate a gap value which is around 7.5 meV 
and 2.7 meV, respectively. Considering the well fixed by 
PCS and tunneling fact |l9| that a gap of about 7 mV 
corresponds to the 2-D a bands and a large gap is seen 
for the direction of the boron planes, we suggest that the 
spectra in Fig. 1 are measured presumably along the de- 
fined direction of the boron sheets. Additional support 
that we have really observed a large gap in this case is 
provided by the magnetic field data. A field about 1 T 
only slightly modified the gap structure in dV/d/ (Fig. 1, 
inset), while according to Ref. , this field is sufficient 
to suppress the small gap. In this connection maxima 
at about 60 mV in the spectra in Fig. 1 are naturally to 
connect with the broadened E2 S phonon modes, which 
predominate in or F(lu) for the a band §. 

Unlike the spectrum in Fig. 1, the spectrum in Fig. 2 
with a small gap about 2.7 meV is attributed to the 3-D n 
bands and is governed by the directions out of the boron 
planes. Correspondingly, the shallow maxima close to 
30 and 50 mV reflect bulk (isotropic) phonons, and their 
correlation with the first two maxima in the phonon DOS 
H or a 2 F(u>) M, H, f?J is evident. We should also empha- 
size that the intensity of the spectrum in Fig. 1 is about 
two times higher than that in Fig. 2, which, according to 
Eq. (2), results in about 5 times larger value for a 2 F{ui) 
due to the difference in the PC resistance Rq. This is 
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FIG. 2: PC spectra for the MgB 2 -Cu contact (R = 1.5 fi, 
Vi(0)=3.2mV) at T = 4.2 K and 5=8 T (solid curves). The 
dashed curves are recorded at zero field and 40 K. Couple 
curves are shown for the both polarities of the applied voltage. 
Dashed straight lines show position of a kink and saturation 
in the PC spectra. The bottom curve is the low energy part 
of the isotropic a 2 F(uj) [fij smoothed by the experimental res- 
olution about 6mV. Inset: the dV /dJ(V) curves for the same 
contact. 
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FIG. 3: PC spectra for MgB 2 -Cu contact (R = 7.5 0, 
Vi(0)=2.5mV) at T = 4.2K in a magnetic field 8T (dashed 
curve) and 9T (solid curve). Curves are almost the same for 
the both polarities of the applied voltage. Inset: dV /dI(V) 
curves for the same contact. 

also in line with an expected larger strength of EPI in 
the boron planes. 

Take a note of the spectra in Figs.l and 2 measured 
above T c around 40 K. Both of them exhibit a buckle 
(kink) at the position of the maximum in the correspond- 
ing low temperature spectra as well as the saturation 
above 80 mV. This is an additional confirmation of the 
phonon caused structure in the curves shown. 

Fig. 3 depicts an extra example of the PC spectrum 
with a smooth maximum around 75 mV which also corre- 
lates with the dominating maximum in many theoretical 
a 2 F{u>) as well as with the Raman data (see Fig. 4(b)). 
Note that though dV /dl is deviated from the BTK shape 
because of a minimum at V—0 and over-gap maxima, 
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FIG. 4: (a) PC spectrum from Fig.l transformed into 
R^dR/dV <x V 2 /V? in comparison with the calculated one 
according to Eq.(3) for the Lorentz-shaped peak at 60 mV 
and a width of 2 mV. The bottom curves show a 2 F(w)|| for 
harmonic and anharmonic phonons. (b) The PC spectrum 
from Fig. 3 in comparison with the calculated one according 
to Eq. (3) for the Lorentz-shaped peak at 70 mV and a width 
of 15 mV. The bottom curve is depolarized Raman spectra of 
MgB 2 @ at 21 K. 

the so-called "gap" minima are observed here at about 
|V| ~10mV. Therefore, we attribute these spectra to the 
rj band. 

Further comparison of our data with the Raman spec- 
trum shows that the Raman maximum (Fig. 4(b)) is nar- 
rower and shifted to a higher energy with respect to that 
in the PC spectrum in Fig. 4(a). Note that the Raman 
spectroscopy studies are restricted to the Brilloum zone 
center, and the Raman maxima, in general, should not 
correspond, e. g., to the maxima in the phonon DOS 
or especially in a.p C F(uj), which reflects phonons with a 
large (maximal) momentum. 

According to the inelastic X-ray scattering data E^] 
an anomalous broadening of the F,2 g modes is found to 
occur along the T-A direction with a width SE =10 - 
17meV and phonon dispersion between 60 and 70meV. 
This energy range hits the maximum in the spectra in 
Fig.l, however, it is fairly broad. Using the width of 
the Fi2 g modes along the T-A direction we estimate the 
phonon lifetime as t a fr/SE ~ (4 — 7) ■ 10~ 14 s. Taking 
into account the upper value of the sound velocity s < 
10 4 m/s |3(| this leads to the phonon mean free path of 
the order of the lattice constant. Therefore, the energized 
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electron generated nonequilibrium phonons accumulate 
in the contact region, giving rise to (i) a large background 
signal 1 25 1 and (ii) a thermalization of phonons directly 
in the PC and a transition to a thermal regime ]25[|. In 
the latter case the phonon features in the PC spectra are 
considerably smeared out. 

The PC EPI spectrum in the thermal regime supposing 
that the phonon part of the resistivity is small compared 
with the residual resistivity is given by Kulik p5| : 



*° av {v) 



c 



duj 



a 2 F(uj)S(eV/huj), (3) 



where C is the constant and Six) represents a smeared 
around x=l step with a shallow maximum at x=1.09. 
We have calculated the PC spectrum by (3) for a 2 F{u) 
where the E2 ff mode was simulated by the Lorentz-shaped 
peak. Comparison with the experimental curves is shown 
in Fig. 4. The results describe fairly well the spectrum 
of the contact from Fig. 3 and shows that the spectrum 
in Fig. 1 has a more distinct maximum, as expected in 
the thermal regime. The position of this maximum cor- 
responds to the main peak in a 2 F{u>) calculated for the 
harmonic phonon model This correlates with the 
statement by Shukla et al. jl2| according to which the 
anharmonic contribution to the phonon linewidth is much 
smaller compare to the electron-phonon one. 

Turn to the EPI inelastic contribution in the PC spec- 
trum, which is about 10% of the total PC resistance 
(see Fig. 1, the inset), it is a few times higher than for 
the MgB2 thin films pPf. However, it is of the same or- 
der as that for the non-superconducting transition metal 
diborides |3lJ|, where the parameter EPI is rather low 
A < 0.1. Conventional explanation considers a deviation 
from the ballistic transport, leading to a decrease in the 
intensity of the PC spectra |25| However, due to the 
strong variation of EPI in MgB2 on the Fermi surface 
and the fact that the large EPI is confined to a small 
volume in the k-space along the T-A direction || the 
PC EPI function ap C F(u) can be vastly different com- 
pared with the thermodynamic a 2 F(uj) one. Therefore, 
to elucidate the details of EPI in MgB2 , the calculation 
of ctp C F(uj) with the mentioned iGfactor accounting for 
strong anisotropy of EPI on the Fermi surface and the 
phonon lifetime effects is very desirable. 

Conclusion. We have measured anisotropic PC spec- 
tra in the single crystal of MgB2. By virtue of the strong 
EPI and the short phonon lifetime a very likely close-to- 
thermal-regime state develops for the most PC's, which 
hinders receiving detailed information about the spec- 
tral EPI function. However, certain of the PC's display 
a spectrum, which corresponds to EPI in the a band, 
and exhibits damped broad maximum just above 60mcV 
caused by the E2 9 phonon modes. The ir band spec- 
trum reveals only lower lying phonon modes. Accord- 
ingly, for the first time the experimental manifestation of 
the anisotropic interaction of electrons with E2 9 phonons 



is presented. Our observation is also in consistent with 
the harmonic model of EPI in MgB2. 
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